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' We describe a quantum computer emulator for a generic, general purpose quantum com- 

puter. This emulator consists of a simulator of the physical realization of the quantum 
0^ ' computer and a graphical user interface to program and control the simulator. We illus- 

trate the use of the quantum computer emulator through various implementations of the 
Deutsch-Jozsa and Grover's database search algorithm. 
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1. Introduction 

I Recent progress in the field of quantum information processing has opened new 

prospects to use quantum mechanical phenomena for processing information. The 
operation of elementary quantum logic gates using ion traps, cavity QED, and NMR 
technology has been demonstrated. A primitive Quantum Computer (QC) [1 — 4] 
and secure quantum cryptographic systems have been build [5 — 7] . Recent theoret- 
ical work has shown that a QC has the potential of solving certain computationally 
hard problems such as factoring integers and searching databases much faster than 
a conventional computer [8 — 13]. 

The fact that a QC might be more powerful than an ordinary computer is based 
on the notion that a quantum system can be in any superposition of states and that 
interference of these states allows exponentially many computations to be done 
in parallel [14]. This intrinsic parallelism might be used to solve other difficult 
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problems as well, such as for example the calculation of the physical properties 
of quantum many-body systems [15 — 18]. In fact, part of Feynman's original 
motivation to consider QC's was that they might be used as a vehicle to perform 
exact simulations of quantum mechanical phenomena [19]. 

Just as simulation is an integral part of the design process of each new gener- 
ation of microprocessors, software to emulate the physical model representing the 
hardware implementation of a quantum processor may prove essential. In contrast 
to conventional digital circuits (which may be build using vacuum tubes, relays, 
CMOS etc.) where the internal working of each basic unit is irrelevant for the 
logical operation of the whole machine (but extremely relevant for the speed of 
operation and the cost of the machine of course) , in a QC the internal quantum dy- 
namics of each elementary constituent is a key ingredient of the QC itself. Therefore 
it is essential to incorporate into a simulation model, the physics of the elementary 
units that make up the QC. 

Theoretical work on quantum computation usually assumes the existence of 
units that perform highly idealized unitary operations. However, in practice these 
operations are difficult to realize: Disregarding decoherence, a hardware implemen- 
tation of a QC will perform unitary operations that are more complicated than those 
considered in most theoretical work. Therefore it is important to have theoretical 
tools to validate designs of physically realizable quantum processors. 

This paper describes a Quantum Computer Emulator (QCE) to emulate var- 
ious hardware designs of QC's. The QCE simulates the physical processes that 
govern the operation of the hardware quantum processor, strictly according to 
the laws of quantum mechanics. The QCE also provides an environment to de- 
bug and execute quantum algorithms (QA's) under realistic experimental condi- 
tions. This is illustrated for several implementations of the Deutsch-Jozsa [20, 21] 
and Grover's database search algorithm [12, 13] on QC's using ideal and more re- 
alistic units, such as those used in the 2-qubit NMR QC [3,4]. Elsewhere [18] 
we present results of a QA to compute the thermodynamic properties of quan- 
tum many-body systems obtained on a 21-qubit hard-coded version of the QCE. 
The QCE software runs in a W98/NT4 environment and may be dowloaded from 
http://rugth30.phys.rug.nl/compphys/qce.htm . 

2. QCE: Quantum Computer Emulator 

Generically, hardware QC's are modeled in terms of quantum spins (qubits) that 
evolve in time according to the time-dependent Schrodinger equation (TDSE) 

i§-^\m) = H{t)\m) , (1) 

in units such that h = 1 and where 
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|^W> = aa,i,...,i;t)U,i,...,i) + a(T,i,...,i;t)|T,i,...,i) + ... 

+ a(T,T,...,T;^)IT,T,...,T) , (2) 

describes the state of the whole QC at time t. The complex coefficients a(J,, J,, . . . , J, 
; t), . . . , a(t, t, • • • , T; completely specify the state of the quantum system. The 
time-dependent Hamiltonian H{t) takes the form [22] 

L 

j,k=l a=x,y,z 
L 

(^J,«,oW + /ij,a,lWsin(/i,at + ¥'j,a))5'° , (3) 

j=\ a=x,y,z 

where the first sum runs over all pairs P of spins (qubits), 5" denotes the a-th 
component of the spin-1/2 operator representing the j-th qubit, Jj,k,ait) determines 
the strength of the interaction between the qubits labeled j and k, hj^a,oit) and 
hj^ce,i{t) are the static (magnetic) and periodic (RF) field acting on the j-th spin 
respectively. The frequency and phase of the periodic field are denoted by fj^a and 
(Pj,a- The number of qubits is L and the dimension of the Hilbert space D — 2^. 

Hamiltonian (3) is sufficiently general to capture the salient features of most 
physical models of QC's. Interactions between qubits that involve different spin 
components have been left out in (3) because we are not aware of a candidate 
technology of QC where these would be important. Incorporating these interactions 
requires some trivial additions to the QCE program. 

A QA for QC model (3) consists of a sequence of elementary operations which 
we will call micro instructions (Mi's) in the sequel. They are not exactly playing 
the same role as Mi's do in digital processors, they merely represent the smallest 
units of operation the quantum processor can carry out. The action of a Ml on the 
state of the quantum processor is defined by specifying how long it acts (i.e. 
the time interval it is active), and the values of all the J's and /I's appearing in 
(3). The J's and /I's are fixed during the operation of the ML A Ml transforms the 
input state \^{t)) into the output state |\l/(t + r)) where r denotes the time interval 
during which the MI is active. During this time interval the only time-dependence 
of H{t) is through the sinusoidal modulation of the fields on the spins. 

Procedures to construct unconditionally stable, accurate and efficient algorithms 
to solve the TDSE of a wide variety of continuum and lattice models have been re- 
viewed elsewhere [23 — 26] . A detailed account of the application of this approach to 
two-dimensional quantum spin models can be found in [27] . Here we limit ourselves 
to a discussion of the basic steps in the construction of an algorithm to solve the 
TDSE for a arbitrary model of the type (3). According to (2) the time evolution 
of the QC, i.e. the solution of TDSE (1), is determined by the unitary transfor- 
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mation U{t + T,t) = exp^.(— z ^ H{u)du), where exp^. denotes the time-ordered 
exponential function. Using the semi-group property of U{t + r, t) we can write 

U{t + T,t) = U{t + m6,t + {m-l)6)---U{t + 26,t + 6)U{t + 6,t) , (4) 

where r = m6 {m > 1) . In general the first step is to replace each U{t + {n + l)d,t + 
nd) by a symmetrized Suzuki product-formula approximation [23, 29]. For the case 
at hand a convenient choice is (other decompositions [27, 30] work equally well but 
are somewhat less efficient for our purposes): 

U{t+ {n+l)6,t + n6) ^U{t+ {n+l)6,t + n6) , (5a) 

U^t +{n+ 1)6, t + nS)^ ^-^SH^{t+{n+l/2)S)/2^-^SH,it+in+l/2)5)/2 

X ^-iSH^{t+{n+l/2)S)^-iSHy{t+{n+l/2)S)/2 
^ ^-i5H,{t+{n+l/2)S)/2 ^ ^g^^ 

where 

L 

Ha{t) = — ^2 Ji,3,aS'jS^ 
j,k=l 
L 

-^{hj,a,o + hj^a,ism{fj^at + (fj,a)) Sj ; a^x,y,z . (6) 

Note that in (6) we have omitted the time dependence of the J's and the /I's to 
emphasize that these parameters are fixed during the execution of a particular MI. 

Evidently U{t + t, t) is unitary by construction, implying that the algorithm 
to solve the TDSE is unconditionally stable [23]. It can be shown that \\U{t + 
T,t) — U{t + T,t)\\ < c6^, implying that the algorithm is correct to second order 
in the time-step 6 [23] . If necessary, U{t + t, t) can be used as a building block to 
construct higher-order algorithms [31 — 33]. In practice it is easy to find reasonable 
values of m such that the results obtained no longer depend on m (and 5). Then, 
for all practical purposes, these results are indistinghuisable from the exact solution 
of the TDSE (1). 

As already indicated above, as basis states we take the direct product 

of the eigenvectors of the Sj (i.e. spin-up | t)j and spin-down | In this basis 

^-iSH^{t+(n+i/2)S)/2 Qiianges the input state by altering the phase of each of the basis 
vectors. As H;g is a sum of pair interactions it is trivial to rewrite this operation 
as a direct product of 4x4 diagonal matrices (containing the interaction-controlled 
phase shifts) and 4x4 unit matrices. Hence the computation of exp{—iSHz{t + {n + 
l/2)5)/2)\^) has been reduced to the multiplication of two vectors, element-by- 
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element. The QCE carries out 0{P2^) operations to perform this calculation but a 
real QC operates on all qubits simultaneously and would therefore only need 0{P) 
operations. 

Still working in the same representation, the action of e-*<5''^a(*+('^+i/2)'5)/2 ^^^^ 
written in a similar manner but the matrices that contain the interaction-controlled 
phase-shift have to be replaced by non-diagonal matrices. Although this does not 
present a real problem it is more efficient and systematic to proceed as follows. Let 
us denote by X (y) the rotation by 7r/2 of all spins about the a;(y)-axis. As 

^-i5Hy{t+{n+l/2)S)/2 ^-iSHy{t+{n+l/2)S) /2 

it is clear that the action of e-'^'^SHyi^+in+^/'^)S)/2 -^^ computed by applying to 
each qubit, the inverse of X followed by an interaction-controlled phase-shift and 
X. The prime in (7) indicates that /ii,z,07 and fi^z in Hz{t + {n + 1/2)5) 

have to be replaced by Jij^y, hi^y^, and fi^y respectively. A similar procedure 

is used to compute the action of e~*''^^(*+'^"'+^/^)''): We only have to replace X 
by y. The operation counts for g-^-^-^- V2)5) (or e-^^Hy{t+{n+i/2)S)/2^ ^^^^ 

0{{P + 2)2^) and 0{P + 2) for the QCE and QC respectively. On a QC the total 
operation count per time-step is 0{3P + 4). 

The operation count of the algorithm described above (and variations of it, see 
e.g. [27]) increases linearly with the dimension D of the Hilbert space, and cannot be 
improved in that sense (although there is may be room for reducing the prefactor by 
more clever programming) . On the other hand one might be tempted to think that 
for small D the cost of an exact diagonalization of the time-dependent Hamiltonian 
(3) at each time-step r may be compensated for by absence of intermediate time- 
steps S. For the problem at hand this is unlikely to be the case: The sinusoidal terms 
in (3) require the use of a time-step that is much smaller than the time-step that 
guarantuees a high accuracy of the Suzuki product formula. Therefore in practice 
for time-dependent Hamiltonian (3), 5 = r so that disregarding implementation 
issues, the running time of the algorithm is very close to the theoretical limit. 

3. Graphical User Interface 

The QCE consists of a QC simulator, described above, and a graphical user in- 
terface (GUI) that controls the former. The GUI considerably simplifies the task of 
specifying the Mi's (i.e. to model the hardware) and to execute quantum programs 
(QP's). The QCE runs in a Windows 98/NT environment. Using the GUI is very 
much like working with any other standard MS-Windows application. The maxi- 
mum number of qubits in the version of the QCE that is available for distribution is 
limited to eight. The QCE is distributed as a self-installing executable, containing 
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the program, documentation, and all the QP's discussed in this paper. These QP's 
also illustrate the use of the GUI itself. 

Some of the salient features of the GUI of the QCE are shown in figs. 1-4. The 
main window contains a window that shows the set of Mi's that is currently active 
and several other windows (limited to 10) that contain QPs. Help on a button 
appears when the mouse moves over the button, a standard Windows feature. 

Writing a QA on the QCE from scratch is a two-step process. First one has to 
specify the Mi's, taking into account the particular physical realization of the QC 
that one wants to emulate. The " MI" window offers all necessary tools to edit (see 
Fig. 2) and manipulate (groups of) Mi's. The second step, writing a QP, consists of 
dragging and dropping Mi's onto a "QP" window. 

Each MI set has two reserved Mi's: A break point (allowing the QP to pause 
at a specified point) and a Ml to initialize the QC. Normally the latter is the first 
instruction in a QP. Each QP window has a few buttons to control the operation 
of the QC. 

The results of executing a QP appear in color-coded form at the bottom of the 
corresponding program window. For each qubit the expectation value of the three 
spin components are shown: = 1/2 — (S'j^) {a = x,y,z) green corresponds to 
0, red to 1. Usually only one row of values (the 2;-component) will be of interest. 
Optionally the QCE will generate text files with the numerical results for further 
processing. 

The QCE supports the use of QP's as Mi's (see Figs. 3, 4). QP's can be added 
to a particular Ml set through the button labeled "QP" . During execution, a QP 
that is called from another QP will call either another QP or a genuine Ml from 
the currently loaded set of Mi's. The QCE will skip all initialization Mi's except 
for the first one. This facilitates the testing of QP that are used as sub-QP's. A 
QP calling a MI that cannot be found in the current MI set will generate an error 
message and stop. 

4. Applications 

Our aim is to illustrate how to use the QCE to simulate the QC implemented 
using NMR techniques [1,2,3,4]. A classical coin has been used to decide which 
of the two realizations (i.e. [1,2] or [3,4]) to take as an example. In the NMR 
experiments [3, 4] the two nuclear spins of the (-"^H and -"^^C atoms in a carbon-13 
labeled chloroform) molecule are placed in a strong static magnetic field in the +z 
direction. In the absence of interactions with other degrees of freedom this spin-1/2 
system can be modeled by the hamiltonian 

H = —Ji,2,zSfS2 — /ii,z,o5'i — /i2,z,o5'| , (8) 
where /ii,2,o/27r !v 500MHz, /i2,z,o/27r 125MHz, and Ji,2,^/27r !^ -215Hz [3]. It 
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is amusing to note that the most simple spin-1/2 system, i.e. the Ising model, can 
be used for quantum computing [34 — 38] 

In the chloroform molecule the antiferromagnetic interaction between the spins 
is much weaker than the coupling to the external field and (8) is a diagonal matrix 
with respect to the basis states chosen, the ground state of (8) is the state with the 
two spins up. Following [3] we denote this state by |00) = |0) (g) |0) = | IT), i-e. the 
state with spin up corresponds to a qubit |0). A state of the A/"-qubit QC will be 
denoted by \x1X2 . . .xn) — \xi) \x2) . . .\xn)- 

It is expedient to write the TDSE for this problem in frames of reference rotating 
with the nuclear spin. Formally this is accomplished by substituting in (1) 

= e**('*i'^'°^i+''2'^'°^2)|*(t)) , (9) 

so that in the absence of RF-fields the time evolution of \E'(t) is governed by the 
hamiltonian H = —Ji,2,zSfS2- 

This transformation removes from the sequence of elementary operations, phase 
factors that are irrelevant for the value of the qubits. Indeed, as the expectation 
value of a qubit is related to the expectation value of the z component of the spin: 

Qj = Q'^=^-{m\S!\m) , (lOa) 

and 

=(vi/(t)|e-**('*i.^.o'^?+'*2,^,oS|)5'|git(hi,^,oSj+/i2,^,oS|)|v[f(^^)^ 
={^{t)\S^\^{t)) . (106) 

^From (10) it is clear that transformation (9) has no net effect. This is not the case 
for the expectation values of the x or y component of the spins: The phase factors 
induce an oscillatory behavior, reflecting the fact that the spins are rotating about 
the 2;-axis (see (^-17) for an example). In the following it is implicitly assumed that 
the basis states of the spins refer to states in the corresponding rotating frame. 

We now discuss the implementation on the QCE of two QA's that have been 
tested on an NMR QC [1, 2, 3, 4]. 
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Table 1. Input and output values of constant {fi{x), f2{x)) and balanced functions (fsix), f4,{x)) 
of one input bit x. 



X 


/i(^) 


f2{x) 




f4ix) 








1 





1 


1 





1 


1 






4.1 Deutsch-Jozsa algorithm 

This QA [20, 39] and its refinement [21] provide illustrative examples of how the 
intrinsic parallelism of a QC can be exploited to solve certain decision problems. 

Consider a function / = f{xi,X2,---,XN) = 0,1 that transforms the N bits 
{xn = 0, 1} to one output bit. There are three classes of functions /: Constant 
functions, returning or 1 independent of the input {xn}, balanced functions that 
give / = for exactly half of the 2^ possible inputs and / = 1 for the remaining 
inputs, and other functions that do not belong to one of the two other classes. Some 
examples of constant and balanced functions are given in tables 1 and 2. 

The Deutsch-Jozsa (D-J) algorithm allows a QC to decide whether a function 
is constant or balanced, given the additional piece of information that functions of 
the type other will not be fed into the QC. For a function of one input variable the 
D-J problem is equivalent to the problem of deciding if a coin is fair (has head and 
tail) or fake (e.g. two heads). In the case of the coin we would have to look at both 
sides to see if it is fair. A QC can make a decision by looking only once (at the two 
sides simultaneously). 



Table 2. Input and output values of constant {fi{x), f2{x)) and balanced functions (f^ix), f4{x), 
fsix)) of three input bits x — {x^, X2, x^}. Note that f4{x) only depends on X2 and is therefore 
rather trivially balanced. 



Xi 


X2 


X3 


fl{xi,X2,X3) f2{xi 




fsixi, X2, 


X3) f4{xi,X2,X3) 


f5{xi,X2,X3) 














1 











1 











1 








1 





1 








1 





1 


1 


1 


1 








1 


1 


1 











1 





1 








1 


1 





1 





1 


1 











1 


1 





1 


1 


1 





1 


1 


1 





1 


1 


1 


1 
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Table 3. Results of letting the sequences F^, . . . , (see (12)) transform the four basis states. 
Inspection of the outputs demonstrated that these sequences implement the constant or balanced 
functions of Table 1. 



X 


1*) 


Fll*) 




F3I*) 


F4I*) 





|00) 


-|00) 


z|01) 


g-z7r/4|00) 


_e+W4|oi) 


1 


|10) 


-|10) 




-e-*''/4|ll) 


e+W4 10) 


2 


101) 


-101) 




g-i7r/4|0l) 


_g + i7r/4 QQ^ 


3 


111) 


-111) 


z|10) 


_e-W4|io) 





In the NMR experiment the two qubits of the QC (i.e. the two nuclear spins 
of the chloroform molecule) are used during the execution of the D-J algorithm 
although in principle only one qubit would do [21]. However our aim is to simulate 
the NMR-QC experiment and therefore we will closely follow Ref. [3]. Accordingly 
the first qubit is considered as the input variable, the other one serves as work 
space. 

Before the actual calculation starts the QC has to be initialized. This amounts 
to setting each of the two qubits to |0). On the QCE this is accomplished by the MI 
"Initialize", a reserved MI name in the QCE (see above). The first step in the D-J 
algorithm is to prepare the QC by putting the first qubit in the state (|0) + |1))/ ^/2 
and the second one in (|0) — |l))/\/2 [3]. This can be done by performing two 
rotations: 

Prepare <^ Y2Y1 , (11) 

where Yj represents the operation of rotating clock- wise the spin j by 7r/2 along the 
y axis, and Yj its inverse (see Appendix A). In this paper we adopt the convention 
that all expressions like (11) have to be read from right to left. 

The next step is to compute the function f{x). Following [3] the two constant 
and two balanced functions listed in Table 1 can be implemented by the sequences 



fi{x) ^Fi= X2X2l{n/2)X2X2l{Ti/2) , (12a) 

/2(x)^F2 = /(7r/2)X2X2/(7r/2) , (126) 

h{x)^F^ = Y^X^Y^X2Y2l{'K)Y2 , (12c) 

/4(a;)^F4 = FiXiYiX2F2/(7r)F2 , (12d) 



where Xj denotes the clock-wise rotation of spin j by 7r/2 along the x axis, Xj 
the inverse operation and /(a) = e"*"*^! represents the time evolution due to H 
itself. In Table 3 we show the result of letting the sequences (12) act on the basis 
states. It is clear that they have the desired properties. Note that prefactors have 
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Table 4. Specification of the micro instructions implementing the two-qubit NMR QC on the QCE. 
Frequencies have been rescaled such that hj^^^o = 1 corresponds to 500MHz. The execution time 
of each micro instruction is given by the second row {t/2tt). The inverse of e.g. is found by 
reversing the sign of i. Note that the QCE is constructed such that a rotation about the x{y) 
axis requires a RF-pulse along the y{x) direction (see Appendix A). 

















Parameter 






Yi 


Y2 


/(7r/2) 


/(tt) 


T/27r 


10 


40 


10 


40 


25 X 10^ 


50 X 10^ 


Jl,2,z 


-10~^ 


-10~^ 


-10~^ 


-10~^ 


-10~^ 


-10~^ 


hi,x,o 




















h2,x,0 








































h2,y,0 




















hi,z,o 


1 


1 


1 


1 


1 


1 


h2,z,0 


0.25 


0.25 


0.25 


0.25 


0.25 


0.25 


hi,x,i 








0.05 


-0.05 








h2,x,l 








0.0125 


-0.0125 








fl,x 








1 


0.25 








f2,x 








1 


0.25 








^l,x 




















f2,x 






















-0.05 


0.05 














h2,y,l 


-0.0125 


0.0125 














h,y 


1 


U.25 


U 











h,y 


1 


0.25 


































'P'2,y 




















hi,z,i 




















h2,z,l 




















fl,z 




















f2,z 




















^l,z 




















V2,z 





















no physical relevance (they drop out when we compute expectation values) and 
that Fi is a rather complicated version of the identity operation. 

Finally there is a read-out operation which corresponds to in the inverse of the 
"Prepare" : 
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Table 5. Final state of the QC after running the D-J algorithm for the case of the ideal QC (Qi, Q21 
see table 5) and the NMR-QC (Qi, O2) see table 4). The results {Qi, O2) have been obtained by 
modifying the NMR Mi's such that the RF-pulses only affect the spin that is in resonance. The 
last two rows show the results of running the refined version [21] of the D-J algorithm. Qi : Ideal 
operations : NMR implementation. 





h{x) 






h{x) 


Qi 


0.000 


0.000 


1.000 


1.000 


Qi 


0.000 


0.000 


0.000 


0.000 


Qi 


0.169 


0.064 


0.867 


0.867 


Q2 


0.999 


1.000 


0.001 


0.002 


Qi 


0.000 


0.000 


0.998 


0.998 


Q2 


1.000 


1.000 


0.001 


0.001 


Qi 


0.000 


0.000 


1.000 


1.000 


Qi 


0.000 


0.000 


0.995 


0.996 



Table 6. Specification of the micro instructions implementing a mathematically perfect two-qubit 

QC on the QCE. The execution time of each micro instruction is given by the second row {t/2'k). 
The inverse of e.g. is found by reversing the sign of /ii,a:,o- Model parameters omitted are zero 
for all micro instructions. 













/(7r/2) 


/(tt) 


Parameter 


Xi 


X2 


Yi 


Y2 


T/27r 


0.25 


0.25 


0.25 


0.25 


25 X 10^ 


50 X 10^ 


-^1,2,2 














-10-6 


-10-6 


^l,x,0 


+1 

















h2,x,0 





-1 














hl,y,0 








+1 











h2,y,0 











-1 









Readout <^ Y1Y2 . (13) 

Note that there is some flexibihty in the choice of these sequences. For instance to 
"Prepare" we could have used Walsh-Hadamard (WH) transformations W1W2 as 
well. 

Upto this point the D-J algorithm has been written as a sequence of unitary 
operations that perform specific tasks. Now we consider two different implemen- 
tations of these unitary transformations: The first one will be physical, i.e. we 
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will use the QCE simulate the NMR-QC experiment itself. The second will be 
"computer-science" like, i.e. we will use highly idealized, non-realizable rotations. 

NMR uses radiofrequency electromagnetic pulses to rotate the spins [40, 41]. By 
tuning the frequency of the RF-field to the precession frequency of a particular 
spin, the power of the applied pulse (= intensity times duration) controls how 
much the spin will rotate. The axis of the rotation is determined by the direction 
of the applied RF-field (see [40, 41] or Appendix A). A possible choice of the model 
parameters, corresponding to the actual experimental values of these parameters, is 
given in Table 4. For simplicity all frequencies have been normalized with respect 
to the largest one (i.e. 500MHz in the experiments [3,4]) . Also note that it is 
convenient to express execution times in units of 27r, the default setting in the QCE. 

The results of running the QCE with the Mi's simulating the NMR experiment 
are summarized in Fig.l. The first qubit {Qj = 1/2 — {Sj)) unambigously tells 
us that the functions fi{x) and f2{x) are constant and that fsix) and /^{x) arc 
balanced. Clearly the QCE qualitatively reproduces the experimental results. In 
the D-J algorithm the final state of the second qubit is irrelevant. In the final state 
the numerical value of qubit 1, is only approximately zero or one (see Table 5). 
This is a direct consequence of the fact that we are simulating a genuine physical 
system. 

In an NMR experiment, application of a RF-pulsc affects all spins in the sample. 
Although the response of a spin to the RF-field will only be large when this spin 
is at resonance, the state of the spins that are not in resonance will also change. 
These unitary transformations not necessarily commute with the sequence of uni- 
tary transformations that follow and may therefore affect the final outcome of the 
whole computation. Furthermore the use of a time-dependent external field to ro- 
tate spins is only an approximation to the simple rotations envisaged in theoretical 
work (see Appendix A). This definitely has an effect on the expectation values of 
the spin operators. 

With the QCE it is very easy to make a detailed comparison between physical 
and idealized implementations of QC's: We simply replace the set of Mi's ("NMR") 
by another one ("Ideal", or "NMR-Ideal") and re-run the QP's by simply clicking 
on the execute buttons. The model parameters we have chosen to implement the 
"ideal" operations are listed in Table 6. The set "NMR-Ideal" is a copy of "NMR" 
(see Table 4) except that the RF-pulses only affect the spin that is in resonance, 
i.e. all operations on qubit j have h2,x,j = h2,yj = 0. 

The results of executing the QA's for the "Ideal" case are shown in Fig. 2. In 
the final state the qubits are exactly |0) or |1), as expected. The state of the second 
qubit not always matches the corresponding state of Fig.l. As mentioned above 
this is due to the approximate nature of the operations used in the NMR case, but 
as the final state of the second qubit is irrelevant for the D-J algorithm there is no 
problem. In Table 5 we collect the numerical values of the qubits as obtained by 
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running the D-J algorithm on the NMR and ideal QC. It is clear in that all cases 
the D-J algorithm gives the correct answer. 



4-2 Collins- Kim- Holton algorithm 

^From the description of the DJ algorithm of one variable ( [3], Fig.l) it is 
evident that the second qubit is redundant because the function call (step T2 in 
[3]) leaves the state of the second qubit, i.e. the work space, untouched. A refined 
version of the D-J algorithm (for an arbitrary number of qubits) that does not 
require a bit for the evaluation of the function is given in [21]. For one variable, 
the QC has to compute the function 

• (14) 
Following [21] this may be accomplished via an /-controlled gate defined by 

Uf\x) = {-iy^-^\x) . (15) 

Accordingly, once choice (there are several) of the set of sequences that implements 
the refined version of the D-J algorithm reads: 

Prepare <^Fi , (16a) 

h{x)^Fr = WrWi , (166) 

f2{x)^F2 = X^Xx , (16c) 

/3(a;)^F3 = yiXiYiYiXiYi , (16d) 

h{x)^F^ = Y^X^Y^Y^X^Y^ , (16e) 

Readout <^yi , (16/) 

The results of running these QP's on the QCE are given in Table 6. It is clear that 
the refined version performs as expected. 



4.-3 Grover's database search algorithm 



On a conventional computer finding a particular entry in an unsorted list of N 
elements requires of the order of N operations. Grover has shown that a QC can 
find the item using only 0{-\fN) attempts [11, 12]. Consider the extremely simple 
case of a database containing four items and functions gj{x), j = 0, . . . , 3 that upon 
query of the database return minus one x = j and plus one if x 7^ j. Assuming 
a uniform probability distribution for the item to be in one of the four locations, 
the average number of queries required by a conventional algorithm is 9/4. With 
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Grover's QA the correct answer can be found in a single query (this result only- 
holds for a database with 4 items). Grover's algorithm for the four-item database 
can be implemented on a two-qubit QC. 

The key ingredient of Grover's algorithm is an operation that replaces each am- 
plitude of the basis states in the superposition by two times the average amplitude 
minus the amplitude itself. This operation is called "inversion about the mean" and 
amplifies the amplitude of the basis state that represents the searched-for item. To 
see how this works it is useful to consider an example. Let us assume that the 
item to search for corresponds to e.g. number 2 ((72(0) = (72(1) = 6*2 (3) = 1 and 
(72(2) = —1). Using the binary representation of integers with the order of the bits 
reversed, the QC is in the state (up to an irrelevant phase factor as usual) 

|*) = 1(|00) + |10)-|01) + |11)) . (17) 

We return to the question of how to prepare this state below. The operator D that 
inverts states like (17) about their mean reads 



/-I 1 1 

2 11-1 
V 1 1 1 

The mean amplitude of (17) is 1/4 and 



1\ 




|00) 


1 




|10) 


1 




|01) 






111) 


find that 





(18) 



we 



D\-<^) =|01) , (19a) 
i.e. the correct answer, and 

d2|*)=^(|00) + |10) + |01) + |11)) , (196) 

^'l*)=-^(|00) + |10)-|01) + |ll)) = -|*) , (19c) 

showing that (in the case of 2 qubits) the correct answer (i.e. the absolute value 
of the amplitude of |10) equal to one) is obtained after 1, 4, 7, ... iterations. In 
general, for more than two qubits, more than one application of D is required to 
get the correct answer. In this sense the 2-qubit case is somewhat special. 

The next task is to express the preparation and query steps in terms of ele- 
mentary rotations. For illustrative purposes we stick to the example used above. 
Initially we set the QC in the state |00), i.e. the state with both spins up [42]. and 
then transform |00) to the linear superposition (17) by a two-step process. First we 
set the QC in the uniform superposition state \U): 
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Prepare ^ \U) = W2Wi\00) = --(|00) + |10) + |01) + |11)) 



where 



(20) 



(21) 



is the WH tranform on qubit j which transforms |0) to z(|0) + |l))/-\/2 (see Appendix 
A). The transformation that corresponds to the apphcation of g2{x) to the uniform 
superposition state is 



/ 1 




V 





1 

-1 








1/ 



|00) 

|10) 
|01) 

111) 



(22) 



This transformation can be implemented by first letting the system evolve in time: 



/(7r)|C/) = e-*'^^^^^^ 



^(|00) + |10) + |01) + |11)) 
1 



= + e+^"^4|10) + e+^'^/^IOl) + e-'"/^|ll)) (23) 

Zi 

For the NMR-QC based on hamiltonian (8) this means letting the system evolve in 
time (without applying pulses) for a time tq = — 7r/Ji,2,z (recall Ji,2,2 < 0). Next 
we apply a sequence of single-spin rotations to change the four phase factors such 
that we get the desired state. The two sequences YXY and YXY (see Appendix 
B) are particulary useful for this purpose. We find 



Y1X1Y1Y2X2Y2 



-(e-^'^/>0) + e+^'^/^llO) + e+^'^/^lOl) + e-^'^/^Ill)) 



^(e-^'^/^O) + e-^^/^|10) + e+^^'^/^lOl) + e-^'^/^Ill)) 

(|00) + |10)-|01) + |11)) . (24) 



2 

g-i7r/4 



Combining (23) and (24) we can construct the sequence Gj that transforms the 
uniform superposition \U) to the state that corresponds to gj{x): 



Fo =YiX,YiY2X2Y2l{n) 
Fi =YiXiYiY2X2Y2l{Tv) 
F2 ^Y^XiYiY2X2Y2l{7i) 

Fs =Y^X^Y^Y2X2Y2l{7T) 



(25a) 
(256) 
(25c) 
(25cZ) 
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The remaining task is to express the operation of inversion about the mean, i.e. the 
matrix D (see (18)), by a sequence of elementary operations. It is not difficult to 
see that D can be written as the product of a WH transform, a conditional phase 
shift P and another WH transform: 



--W1W2 



W1W2 . (26a) 



D =WiW2PWiW2 

/ 1 0\ 
0-100 
0-10 

V -1 / 

The same approach that was used to implement g2{x) also works for the conditional 
phase shift P (= —Fq) and yields 

P = YiX^nY2X2Y2l{n) . (27) 
The complete sequence Uj reads 

Uj = WiW2PWiW2Fj . (28) 

Each sequence Uj can be shortened by observing that in some cases a rotation is 
followed by its inverse. Making use of the alternative representations of the WH 
transform Wi (see (Appendix B)), the sequence for e.g. j = I can be written as 



WiW2Fi = - XiXiYiX2X2Y2YiX^YiY2X2Y2l{'K) 

= - XiYiX2Y2l{7r) . (29) 

The sequences for the other cases can be shortened as well, yielding 

Uo^XiYiX2Y2l{T^)XiYiX2Y2lin) , (30a) 

Ui ^XiYiX2Y2HT^)XiYiX2Y2l{^) , (306) 

U2 =XiYiX2Y2I{tt)XiYiX2Y2I{tt) , (30c) 

Us =XiYiX2Y2l{7r)XiYiX2Y2l{7r) , {30d) 

where in Ui and t/2 a physically irrelevant sign has been dropped. Note that the 
binary representation of x translates into the presence (0) or absence (1) in (30) of 
a bar on the rightmost Xi and X2. 

As before, our aim is to use the QCE to simulate the NMR-QC experiment 
[4]. For the D-J algorithm we already specified the physical parameters for the 
elementary operations and we will make use of the same set of Mi's here. In Figs. 3 
and 4 wc show the QCE after running the four cases go{x), . . . , g3{x) using the NMR 
(Fig. 4) Mi's. The numerical values of the qubits in the final state are given in Table 
7, for the ideal and NMR QC. In both cases the QA performs as it should. In the 
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Table 7. Final state of the QC after running the Grover's database search algorithm for the case 
of the ideal QC (Qi, Q2, see table 5) and the NMR-QC (Qi, Qa, see table 4). 





9o{x) 


9i{x) 


92{x) 


93{x) 


Qi 





1 





1 


Q2 








1 


1 


Qi 


0.028 


0.966 


0.037 


0.955 


Q2 


0.163 


0.171 


0.836 


0.830 



ideal case, the final state of the QC is exactly equal to \x) (binary representation of 
integers). Using RF-pulses instead of ideal transformations to perform 7r/2 rotations 
leads to less certain answers: The final state is no longer a pure basis state but some 
linear superposition of the four basis states. What is beyond doubt though is that 
in all cases the weight of \x) is by far the largest. Hence the QC returns the correct 
answer. 

5. Summary 

We have described the internal operation of QCE, a the software tool for simulat- 
ing hardware realizations of quantum computers. The QCE simulates the physical 
(quantum) processes that govern the operation of the hardware quantum processor 
by solving the time-dependent Schrodinger equation. The use of the QCE has been 
illustrated by several implementations of the Deutsch-Jozsa and Grover's database 
search algorithm, on QC's using ideal and more realistic units, such as those of 
2-qubit NMR-QC's. Currently the QCE is used to study the stability of quantum 
computers in relation to the non-idealness of realizable elementary operations [43] . 
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Appendix A: Spin- 1/2 algebra 

Here we present a collection of standard results on spin- 1/2 systems which are 
used in the paper and are taken from [41]. We begin with some notation. 

The two basis states spanning the Hilbert space of a two-state quantum system 
are usually denoted by 
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It). (J) , 11). (°) . (Al) 

The three components of the spin- 1/2 operator S acting on this Hilbert space are 
defined by 

By convention the represenation (A. 2) is chosen such that | t) and | |) are eigen- 
states of with eigenvalues +h/2 and —h/2 respectively. 

^From (A.2) it is clear that {S'^f = {Syf = {S^f = h'^/A so that 

003(5^) =coK|) (J I) . (A3a) 

and 

sin(^)=|sin(|)5" . {AM) 
The commutation relations between the three spin-components read 

where [A, B] = AB — BA, e^pry is the totally antisymmetric unit tensor ( e^yz = 
^yzx = ^zxy = I5 ^afi^y = —^pa-y = —^ji3a = —£0:7/35 ^007 = 0) and the Summation 
convention is assumed. 

Rotation of the spin about an angle (f around the axis P gives 

S'^{ip,P) = e''^^'/''S'^e-''^^'/^ = 3'' cos ip + e^p^S'' sin if . (A.5) 

Of particular interest to quantum computing are rotations about tt/2 around the x 
and y-axis defined by 

and 

Y^e-''/^^ = ±(^_l J) . (AM) 

The inverse of a rotation Z will be denoted as Z and if more than one spin is 
involved a subscript will be attached. With our convention (f |yS'^y| t) = —1/2 
so that a positive angle corresponds to a rotation in the clock-wise direction. 

Another basic operation is the Walsh-Hadamard transform W which rotates the 
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state I t) into (| + | (up to an irrelevant phase factor), i.e. the uniform su- 

perposition state. In terms of elementary rotations the Walsh-Hadamard transform 
reads 

For example 

'^iT)=w'(j)-7f(;) ■ (^-s) 

We now consider the time evolution of a single spin subject to a constant mag- 
netic field along the 2;-axis and a RF-field along the x-axis, i.e. the elementary 
model of NMR. The TDSE reads 

in^\^t)) = - [HqS^ + HiS^'suiiot] , (y4.9) 

ot 

where = 0)) is the initial state of the two-state system and we have set the phase 
in (3) to zero for notational convenience. Substituting = e**'^0'^'/^|*(t)) 

yields 

in — \-^(t)) = - [{Ho - uo) S"" + Hi S"" sin ujt cos ujot + Hi sin ujt sin ujot] , 

(AlO) 

which upon chosing loq = Hq can be written as 

ih^\^ {t)) = - Hi [S'' sin ojt cos Hot + sin cut sin Hot\\^{t)) . (All) 
At resonance, i.e. co = Ho, we find 

ih^l^it)) = -^[Sy + S'=sin2Hot-Sycos2Hot]\^(t)) . {A.12) 

Assuming that the effects of the higher harmonic terms (i.e. the terms in sin 27^0^ 
and cos 2Hot) are small [41] we obtain 

, (A13) 

which is easily solved to give 

\^{t)) ^(^^^'^'/^^\^{t^Q)) , (A.14) 
so that the overall action of an RF-pulse of duration r can be written as 
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(A15) 



^Prom (A. 15) it follows that application of an RF- pulse of "power" tHi = tt will 
have the effect of rotating the spin by an angle of 7r/2 about the j/-axis. For example 



JrHoSyn inS^ /2h\ 



T) 



TT / 1 0\ . . TT /O 

^^^4 1 +^^^"4U 



V2 

I / giTHo/2 

In this rotated state the expectation values of the spin components are given by 



{A.16) 



(t \e-i^sy/2h^-irHoSyhsx^irHoSyh^insy/2h^ 1^ ^ _ ficosTHo , {A.17a) 



(t \e-i^sy/2h^-irHoSyhgy^irHoSyh^i^sy/2h^ l^ ^ _ fismrHo , (A.176) 



^1 ^^-iirSy /2n^-iTHoSyhgz^iTHoSyh^inSy /2h^ _q 



(A.llc) 



showing that the time of the RF-pulse also affects the projection of the spin on the 
X and y axis. 

It is instructive to derive the TDSE that corresponds to approximation (A. 15). 
Taking the derivative of (A. 15) with respect to t we obtain 



ih^mt)) = -HoS' - ^ [S"" sin Hot + cos Hot] |$(t)) 



(A.18) 



telling us that the approximate solution (A. 15) is the exact solution for an RF 
field rotating in space [41] . The fact that the application of an RF-pulse does not 
exactly correspond to a simple rotation in spin space may well be important for 
applications of NMR techniques to QC's. 

Finally we note that our choice of using a "sinc<;t" instead of "cos cut" [41] to 
couple the spin to the RF- field merely leads to a phase shift. In the former case 
rotating the spin around the a;-axis requires a pulse along the y-axis, whereas in 
the latter the pulse should be applied along the a;-axis [41] . 
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Appendix B: Basic operations 



Below we list a number of identities that are useful to compute by hand the 
action of the sequences appearing above. The convention adopted in this paper is 
that 

|0) = |T) = (J 



|0) = |T) = (J) ; |i) = U) 

A straightforward calculation yields: 

x\o)=^m + i\i)) , x|i) = i= 



'0^ 



x|o)=i=(|o) + i|i)) , x|i) = i=(i|o) + |i)) 

^|0)=^(|0) 

y|o)=i=(|o)-|i)) , F|i) = ^(|o) + |i)) 



F|0)=i=(|0) + |l)) 



Y 
Y 



^(|o)-|i)) 



Y 



;^(|o> + |i)) 



|0) 



Y 



1 



[V2 



:(|0) + |1))] =|1> 



= |0) : 

XX\0)=i\l) , XX\l) = i\0) 

YY\o) = -\i) , yy|i) = |o) 

^>^|0)=^(|0)-|1)) , XY\1) = —^ 

g-i7r/4 



■(|0) + |1)) : 

-\-i7V / 4: — ZTT /4 

^>^|0)=^(|0)-|1)) , XF|1) = ^(|0) + |1)) . 

p+i7r/4 _ p-i7r/4 

.(|0) + |1)) , XF|l) = -^(-|0) + |l)) 

■|o) + |i)) 



XY\0) : 



g+i7r/4 



J— i7r/4 

1^(|0) + |1)) , XY\1) 

YXY\0) =e+'''/'^\0) , = e-^'^/^ll) , 

Fxy|0) =e-^''/'^|0) , YXY\1) =e+'''/'^\l) , 

YXY\0) =e+*^/^|l) , YXY\1) = -e-'^'/^lO) 

YXY\0) =e-'''/^\0) , YXY\1) =e+'''/^\l) , 



(B.l) 
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W =X^Y = YX^ = -X^Y = -YX^ , {B.8a) 

W|0)=-^(|0) + |1)) , W^|1) = -^(|0)-|1)) , (BM) 
^2|o) = _|o) , W^\l) = -\1) . (5.8c) 
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Fig.l. Picture of the Quantum Computer Emulator showing a window with a 
set of micro instructions implementing an NMR quantum computer and windows 
with four Deutsch-Jozsa programs (d-jl, d-j4), one for each function 

fiix)) listed in Table 1. The final state of the QC, i.e. the expectation value 
of the qubits (spin operators), is shown at the bottom of each program window 
(green = |0), red =|1)). The numerical values appear if the cursor moves over the 
qubit area. 
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Fig. 2. Picture of the Quantum Computer Emulator showing a window with a 
set of micro instructions implementing an ideal quantum computer and windows 
with four Deutsch-Jozsa programs (d-jl, d-j4), one for each function {fi{x), 
f^ix)) listed in Table 1. Also shown is a window for editing micro instructions, 
which appears by double-clicking on a micro instruction (x2 in this example) . The 
final state of the quantum computer, i.e. the expectation value of the qubits (spin 
operators), is shown at the bottom of each program window (green = |0), red 
= |1)). In this ideal case the expectation values are either zero or one. 
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Fig. 3. Picture of the Quantum Computer Emulator showing a window with a set 
of micro instructions for an two-qubit NMR quantum computer and windows with 
quantum programs implementing Grover's database search for the four different 
cases goix) (gO), (73 (x) (g3). This example also shows the use of quantum 
programs as micro instructions in other quantum programs. The final state of the 
QC, i.e. two qubits shown at the bottom of each program, gives the location (in 
binary representation) of the item in the database. Note that for the case gl we 
were using single-step mode to execute the program and stopped at fl. 
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Fig. 4. Picture of the Quantum Computer Emulator showing a window with a set 
of micro instructions for the two-qubit NMR quantum computer and windows with 
quantum programs implementing Grover's database search for the four different 
cases (7o(x),... (73(x), using the basic Mi's (grovO, grov3) and calls to other 
quantum programs (gO, g3). The final state of the quantum computer, i.e. 
two qubits shown at the bottom of each program, gives the location (in binary 
representation) of the item in the database. This final state is no longer a pure 
basis state but as the weight of basis state corresponding the location of the item 
is by far the largest the correct answer is easy to infer. 



